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Abstract 

We study zero-range processes which are known to exhibit a conden- 
sation transition, where above a critical density a non-zero fraction of 
all particles accumulates on a single lattice site. This phenomenon has 
been a subject of recent research interest and is well understood in the 
thermodynamic limit. The system shows large finite size effects, and we 
observe a switching between metastable fluid and condensed phases close 
to the critical point, in contrast to the continuous limiting behaviour of 
relevant observables. We describe the leading order finite size effects and 
establish a discontinuity near criticality in a rigorous scaling limit. We 
also characterise the metastable phases using a current matching argu- 
ment and an extension of the fluid phase to supercritical densities. This 
constitutes an interesting example where the thermodynamic limit fails 
to capture essential parts of the dynamics, which are particularly relevant 
in applications with moderate system sizes such as traffic flow or granular 
clustering. 

keywords: zero range process; condensation; metastability; finite size effects; 
large deviations 



1 Introduction 

Zero-range processes are stochastic particle systems with no restriction on the 
number of particles per site and with jump rates that depend only on the oc- 
cupation of the departure site. This simple zero-range interaction leads to a 
product structure of the stationary distributions [37J Q] . These processes have 
been a focus of recent research interest since they can exhibit a condensation 
transition. This is the case for space homogeneous jump rates g(n) that de- 
cay asymptotically with the number of particles n. A prototypical model with 
g(n) = 1 + for n = 1, 2, . . . has been introduced in [12], where conden- 
sation occurs for parameter values 7 G (0, 1), b > or 7 = 1, b > 2. If the 
particle density p exceeds a critical value p c , the system phase separates into a 
homogeneous background with density p c and a condensate, where the excess 
particles accumulate on a single randomly located lattice site. This transition 
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has been established on a rigorous level in a series of papers, in the thermo- 
dynamic limit HU 21 [3J , as well as on a finite system as the total number 
of particles diverges |16j . Dynamic aspects of the transition such as equilibra- 
tion and coarsening [TTJ [3T] and the stationary dynamics of the condensate [T3] 
are well understood heuristically, for the latter first rigorous results have been 
achieved recently [7J [B] . Findings for the zero-range process could be applied 
to understand condensation phenomena in a variety of noncquilibrium systems 
(see [13] and references therein), as well as providing a generic model of do- 
main wall dynamics and a criterion for phase separation using a mapping to 
one-dimensional exclusion systems [27]. The process continues to be of interest, 
recent work on variations of the model includes mechanisms leading to more 
than one condensate [Ml US] j or the effects of memory in the dynamics [33J . 

While most of the results so far consider the thermodynamic limit, finite size 
effects in the model with jump rates g(n) and 7=1 have been investigated in 
[Pol 155] using saddle point methods, and in [5] for a variant of the model with 
a single defect site. In the condensed phase region, finite systems are found 
to exhibit a large overshoot of the stationary current above its value in the 
thermodynamic limit. In this paper we examine this phenomenon in detail for 
all possible values of the parameter 7 € (0, 1]. We find the leading order finite 
size effects that describe the current overshoot by continuing the homogeneous 
(fluid) phase above the critical density and characterizing the condensed phase 
by a current matching argument. For 7 < 1, the main focus of this article, 
finite systems exhibit a metastable switching behaviour between the two phases, 
which is prevalent in Monte Carlo simulations for a wide range of parameters. 
To capture this phenomenon we examine the system in a scaling limit and 
derive a rate function, which exhibits a double well structure, describing the 
distribution over the bulk density (the bulk density serves as an order parameter 
to distinguish the fluid and condensed phases). This way we rigorously establish 
the discontinuous behaviour on the critical scale, even though the bulk density 
is a continuous function in the thermodynamic limit. This is shown to be in 
agreement with recent results on the condensation transition at the critical 
density [3j. Based on the exact scaling limit, we can also predict the lifetime 
of the metastable phases for large finite systems by a heuristic random walk 
argument. 

In general, finite system size can lead to effective long-range interactions 
and non-convexity of thermodynamic potentials, as has been observed for vari- 
ous models (see for example [25] [5] and references therein) . For the zero-range 
process it has been shown that for simple size-dependent jump rates metasta- 
bility effects can be manifested even in the thermodynamic limit [5D] , and large 
crossover effects in these systems have already been observed in [T2]. The 
finite-size behaviour of the zero-range process considered here also exhibits the 
above non-convexity, with the additional feature of a sharp crossover between 
a putative fluid and condensed phase in the non-convex part. This leads to a 
metastable switching behaviour which disappears in the thermodynamic limit, 
intriguingly contradicting the usual expectation of finite systems to behave in a 
smoother fashion than the limiting prediction. The onset of phase coexistence 
at criticality is a classical question of general interest in phase separating sys- 
tems, see for example [9] HOj for the formation of equilibrium droplets in the 
Ising model which also form suddenly on a critical scale. Rigorous results on 
metastability regarding the dynamics of the condensate in zero-range processes 
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have also been a subject of recent research interest |H1 H], and our work pro- 
vides a contribution in that direction and new insight in the mechanisms of 
condensate dynamics on finite systems. This is explained in more detail in the 
discussion. 

A proper understanding of the metastability phenomenon exhibited by zero- 
range processes on finite systems is also of particular importance for recent ap- 
plications with moderate system sizes. Clustering phenomena in granular media 
can be described by zero-range processes (see [3H] [31] and references therein) , 
and metastablc switching between homogeneous and condensed states has been 
observed experimentally [40] [33] . In the spirit of the mapping introduced in [27] 
the zero-range process with jump rates g(n) has also been applied as a simpli- 
fied traffic model [55] [30] , where condensation corresponds to the occurance of 
a traffic jam. A key feature of traffic models is the existence of a broad range 
of densities over which metastability between free flowing and jammed states 
is observed (see for example [H] and references therein). The relevance of this 
study in application serves as a motivation, but the aim of the paper is a general 
understanding of finite size effects and their implications for a generic class of 
zero-range processes, rather than a detailed analysis of particular cases. 

The paper is organized as follows: In Section 2 we introduce the model and 
summarise previous results on stationary distributions and the thermodynamic 
limit. In Section 3 we present relevant observations on finite systems which we 
later study in detail. In Section 4 we give preliminary results and a heuristic 
description of the finite size effects. This motivates our main results presented 
in Section 5, where we rigorously establish metastability in a suitable scaling 
limit. In Section 6 we connect these results to the lifetimes of the metastable 
phases, and end with a short discussion in Section 7. 

2 Stationary Measures and the Thermodynamic 
Limit 

2.1 The Zero Range Process 

We consider a one dimensional lattice of L sites Al = {1, 2,...,L} with periodic 
boundary conditions. Let r\ x £ No = {0, 1,2,.. .} be the number of particles on 
site x € Al. The state of the system is described by r\ = (i] x )xeA L belonging to 
the state space of all particle configurations fit = Nq l . Particles jump on the 
lattice at a rate that depends only on the occupation number of the departure 
site. A particle jumps off site x after an exponential waiting time with rate g{rj x ) 
and moves to a target site y according to the probability distribution p(y — x). 
We assume that p is of finite range, i.e. p(z) = if \z\ > R for some R > 0, 
normalised and irreducible on A^. The main results of the paper focus on the 
jump rates g : No — > K + of the form 



with 7 e (0, 1) and b > 0. A discussion on extending these results to the case 
7 = 1 and b > 3 and other lattice geometries can be found in section [7] These 
jump rates were first introduced by Evans |12j and represent a fairly general 
class of functions of interest for the condensation transition. 
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The infinitesimal generator of the process acting on suitable test functions 
/ is given by 

{Cf){n)= E 9(v*)p(y-x)(f(v x ' v )-f(v)) , (2) 

x.yeA L 

where rj^' y = r\ z — S(z,x) + S(z,y) and 5 is the Kronecker delta [T] EI]. The 
process conserves the total number of particles in the system, so £Il can be 
partitioned into invariant subsets Ql,n = {v ^ E^eA Vx ~ -W} on which 
the zero-range process is a finite state irreducible Markov process. The process 
can also be defined on an infinite lattice under certain constraints, for details 
see [22 E]. 

2.2 Stationary Measures 

The following summarises well known results on stationary measures of the zero- 
range process, for details see [371[T21[T]. The zero-range process with generator 
^ has a family of stationary homogeneous product measures on fi^ which we 
refer to as the grand canonical ensemble. These measures are parameterized 
by a fugacity <fi and are of the form, 

"£M= II v <t>Vlx\ where i> [n] = — l —w(n)cf> n . (3) 

xeA L ^P> 

These exist for all 4> <E [0, </> c ) where 4> c is the radius of convergence of the (single 
site) partition function 

oo 

z(0) = 5>(fc)0 fc . (4) 

fe=0 

The stationary weights w are given by w(0) = 1 and 

n 

w(n) = Y[g(k)-\ n>0. (5) 

fc=i 

In the grand canonical ensemble the expected particle density is a function 
of 4> and is given by, 

oo 

R{4>) :=E^[ m ) =£>^(fc) = Wog*#) , (6) 

fc=0 

which is strictly increasing and i?(0) = 0. The critical density is defined by 
p c = lim R(4>) G (0, oo], and condensation occurs if p c < oo, as explained 

<t>S<f>a 

below. 

The expected jump rate off a site is proportional to the average station- 
ary current or, in case the first moment ^2 z zp(z) vanishes, to the diffusivity. 
Therefore for simplicity, in the rest of this paper, current will refer to the av- 
erage jump rate off a site, which is clearly site independent under a stationary 
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distribution. In the grand canonical ensemble the current is simply given by the 
fugacity tf>, 

34, := E„, [ff(Tte)] = 7-T E «(»W n )^ = ^ . ( 7 ) 
„_ 

which follows directly from the form of the stationary weights u>(n) Eq. (|5j). 

For fixed L and iV the process restricted to FIl,n is ergodic, the correspond- 
ing unique stationary measures belong to the canonical ensemble and are 
given by 

ttl.jv [v] : = v <i> [v\Sl{v) = N] where S L (r]) = ^ rj x . (8) 

x€A L 

These measures are independent of <f> and are given explicitly by 

n L ,N [v] = Z , L jy-N II ^fe^CE^'^) ' 
where the canonical partition function is 

In the canonical ensemble the form of the stationary weights, Eq. ([5]), imply 
that the average current is given by a ratio of partition functions, 

. , Z(L,N- 1) , , 

3l,n ■= E nL N [g( Vx )} = V Z(L (10) 

For the jump rates ([!} the single site weights asymptotically decay as a 
stretched exponential 

w{n) ~ e^T^™ 1 ^ for 76 (0,1). (11) 

Throughout the paper we use '~' to mean asymptotically proportional and '~' 
to mean asymptotically equal. It follows that <j> c = 1 and the critical density 
and the variance are finite, 

p c := R{1) < 00 , a 2 c := [rf x ] - p\ < 00 . 

i? is strictly increasing, so invertible on [0, <fi c ] = [0, 1] and we denote its inverse 

by 

H>(p)=R-\p). (12) 

In this way we can parameterise the grand canonical measures by densities 
pe [0,Pe]- 
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2.3 Thermodynamic Limit 

In this section we summarise known results on the condensation transition in 
the zero-range process. It has been established in [T2] and rigorously in [3TJ [TO] 
as a continuous phase transition in the thermodynamic limit, as particle number 
N and lattice size L tend to infinity such that N/L — > p. It was shown that 
all finite dimensional marginals of the canonical measure ttn,l converge to the 
grand canonical measure with density p if p < p c . If p > p c then there is no 
grand canonical measure with density p and all finite dimensional marginals of 
kl,n converge to the grand canonical measure with density p c . This result on 
the equivalence of ensembles holds in terms of weak convergence (for details see 
|19j ) and can be summarised as 

tl,n vq,(p\ as L — >• oo and N/L p , (13) 

where 

l <Pc if P > Pc 

This implies that for p > p c the excess particles in the system condense on a 
set of vanishing volume fraction. The result has been strengthened in [H [3] 
(with partial results already in [2H]), showing that the condensate typically 
resides on a single lattice site. Denoting the size of the maximum component 
Ml (if) = max r\ x this can be written as 

\ Ml ^{ ° l P ~ Pc . (15) 

L I p — p c h P> Pc 



Here the notation denotes converges in probability with respect to ttl,n, 

for all e > 0. 



~M L -(p-p c ) 



> e 



It will often be useful to express this result also in terms of the behaviour of 
the bulk of the system. We denote the total number of particles outside of the 
maximally occupied site by S L s (rj) = N — M^rf). The background density 
converges as 

^M^U 'Zl% ■ <*> 

A similar result holds for the behaviour of the canonical current in the thermo- 
dynamic limit. With Eq. ([7]) the convergence in (|13p . ([T4"]) implies that 

jh,N^^(p) as L ->• oo and N/L^-p. (17) 

So the current and the background density are both continuous with respect 
to the total system density in the thermodynamic limit. Also both are strictly 
increasing up to p c and constant for p > p c . If N/L — > p < p c the system is 
said to be in the fluid phase region and if N/L — > p > p c the system is in the 
condensed phase region. 
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Density p Density p 

Figure 1: Finite size effects and current overshoot for 7 = 0.5 and 6 = 4. 
Left: The canonical current for various system sizes as a function of the density 
p = N/L are plotted. The dashed black line shows the thermodynamic current 
as a function of the system density p. Right: The overshoot region for L = 1000 
showing the two distinct currents measured from Monte Carlo simulations for 
various densities near the maximum current. 



The thermodynamic entropy is defined by the Legendre transform 

s{p)= sup (log z ((/>)- p log <j>) 
4>e[o,4>a) 

= log*($(p))-plog$(p) , (18) 

where $(p) is given by Eq. (fl"4"]). It has been shown in [2T] that the canonical 
partition function converges to the thermodynamic entropy, i.e. 

- log Z(L,N) -> s(p) asi^oo and N/L^p. (19) 

3 Observations on a Finite Systems 

In this section we present results obtained from exact numerics and Monte Carlo 
simulations in the canonical ensemble. We can calculate the canonical current 
given by Eq. (fTU)) by making use of the following recursion relation for the 
canonical partition functions, 

N 

Z(L,N)=^2w(k)Z(L-l,N-k). (20) 

k=0 

Similarly we can calculate the canonical distribution of the maximum site occu- 
pation Ml ■ To this end we define the cut-off canonical partition function which 
counts configurations for which Ml (17) < m, 

min{m,7V} 

Q(L,N,m)= w(k)Q(L-l,N -k,m) . (21) 

This allows us to calculate 

rw 1 Q(L,N,m)-Q(L,N,m-l) 
n L ,N [Ml = m] = , (22) 
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Figure 2: Finite size effects for the background density for 7 = 0.5, b = 4 and 
L = 1000. Left: The expected value of the background density. The dashed 
black line shows the thermodynamic limit result. Right: Distribution of the 
background density at three system densities shown by corresponding dash- 
dotted lines on the left, calculated exactly using Eq. (|2"2"j) . The position of local 
maxima of the distributions are marked on both plots ( < for p = 0.6, x for 
p = 0.7, o for p = 0.8). The high background density maximum at p = 0.8 
occurs with extremely low probability and is off the scale. 



for all m E N (the case to = is trivial). We often consider the equivalent 
formulation using background densities 



ttl,n [M l = m] = tt l 



n 



1 



v, N - 
where ok, 



L - 1 



This is more useful for illustrations and is more intuitive, since the background 
density characterises all but a single lattice site, while the formulation involving 
Ml is more convenient for computations since it avoids issues with non-integer 
numbers. Therefore keeping both formulations in parallel is the best option for 
a concise presentation of our results. 

On large finite systems we observe significant finite size effects above the 
critical density p c , Fig. [T] shows the typical behaviour of the canonical current. 
Below p c the current is very close to the thermodynamic limit result even on 
relatively small systems (L ~ 100) and the leading order finite size effects can be 
understood immediately from the proof of the thermodynamic limit result |21j . 
However the canonical current and background density significantly overshoot 
their critical values see Fig. [T] and Fig. [5] An overshoot has been observed 
before for systems containing a single site defect [5] for 7 = 1. For 7 < 1 the 
current increases monotonically with density p = N/L, in a way that appears to 
vary only very slightly with system size, up to some size dependent maximum 
current at ptrans(L). For p < ptrans(i) the background density in the system 
is typically very close to p (Fig. For a fixed system size L we associate the 
region p < ptrans(£) with a putative fluid phase. Close to p = ptrans(-k) there 
is an abrupt decrease in the current and background density. For p > ptrnns(L) 
the current and background density are both decreasing in p and tend to their 
respective critical value as p — > 00. We associate the region p > ptrans(^) with 
a putative condensed phase since the increasing density is entirely taken up 
by a large number of particles condensing on a single lattice site. The size of 
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Figure 3: Switching dynamics for 7 = 0.5, 6 = 4 and L = 1000 taken at 
N/L = ptrans(lOOO) = 0.695. Left: (Bottom) Current against time from Monte 
Carlo simulations, calculated by taking the average number of jumps in the 
system over small time windows. The black line indicates the transition between 
the two putative phases. (Top) The location of the maximum does not change 
until the system is fluid. Right: Cumulative tail of the distribution of waiting 
times in the putative fluid phase on a log linear scale. The solid line shows an 
exponential fit. 



the effective fluid overshoot increases with increasing 6 and with decreasing 7, 
and is already very pronounced at 7 = 0. 5, 6 = 4. 

For p = N/L close to ptrans(i) we observe that the canonical distribution 
over background densities has two maxima of similar magnitude (Fig. [2]) and 
Monte Carlo simulations show that the system switches between the two puta- 
tive phases. Fig. [3] shows the typical behaviour of the current close to ptrans{L) 
as a function of time, the putative phases can be clearly distinguished by the 
current. For p w Ptrans(i) in the condensed phase the location of the conden- 
sate does not change, while its position fluctuates heavily in the fluid phase, 
supporting the fact that the particles are distributed homogeneously. 

The empirical distribution of waiting times in the two putative phases is very 
close to an exponential (Fig. fright). This suggests that the switching process is 
approximately Markovian over the range of parameters and jump distributions 
observed, and constitutes a genuine metastability phenomenon. The rate of the 
switching depends on the parameters 6 and 7 as well as the jump distribution 
p(x), which will be discussed in Section [B] in more detail. 



4 Preliminary results 

4.1 Heuristics and Current Matching 

As suggested by the form of the current overshoot in Fig 1, our approach is to 
approximate the putative fluid phase by extending the grand canonical current 
above </> c , which we achieve by means of a cut-off grand canonical measure. Since 
the total number of particles is fixed canonically to N the system can not explore 
states where any single site contains more than TV particles. We therefore expect 
the distribution of particles in the background under the canonical measure will 
always be closer (in any reasonable sense) to a grand canonical measure, with 
some suitably chosen cut-off and fugacitiy chosen to fix the correct density of 
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particles, than it is to the unconditioned distribution. For cut-off m G N the 
cut-off grand canonical measures are defined by single site marginals with 
support on {0, 1, . . . , m}, 

v<S> t m Vix = n] := v$ \r\ x = n\rj x < m] 

w{n)(j) n for n < m (23) 



z m {4>) 

where the normalisation is given by the finite sum 



*mO) =5Z^(fc)0 fc . (24) 

These measures are well defined for all <j> € (0, oo) and the current is given by 

^ r / m , / w(m)(b' m \ 

Also the average density in the cut-off ensemble 

1 m 

i? m (0) := E^ m [ Vx ] = —— J2 kw(k) (26) 

m\<P) fe=Q 

is a strictly increasing function from [0, oo) onto [0,oo) and so we denote its 
inverse 

*m(pbg) := R^ipbg)- (27) 

To a first approximation we estimate the current in the putative fluid phase 
using Eq. fl2S), 

3L.N « *Jv(p) (28) 

under the assumption that (l — w(N)(f> N / zn((/))^ w 1 for values of under 
consideration. This approximation is shown in Fig. |4] by the blue dashed line 
for 7 = 0.5, b = 4 and a system size of L = 1000. The approximation is 
extremely close to the canonical current for p < ptva.ns{L) and for p « Ptrans(^) 
is in good agreement with the empirically measured fluid currents from Monte 
Carlo simulations. We demonstrate in the Section [5] that this estimate can be 
improved by choosing the cut-off more carefully. 

We approximate the current in the putative condensed phase by a current 
matching argument, a similar heuristic argument has been given before in [28j . 
The existence of a stable condensate implies that the average rate of particles 
exiting the condensate must be equal to the average rate of particles entering 
it. Conditioned on the occupation of the condensate Ml = m, the exit rate is 
simply g(m) while the entry rate is well approximated by the stationary current 
in the background. This is assumed to be in a putative fluid phase described 
by ^0. m , which leads to the current matching condition. 

N — 7TL 

$m (pb s ) = 9{ m ) where p hg = — — . (29) 

Li — 1 

The lower branch of solutions to this equation define the condensed current ap- 
proximation (Fig. [4] red dashed line). These are extremely close to the canonical 
current for p > ptrans(^) and for p sa ptrans(^) are in good agreement with the 
empirically measured condensed currents from Monte Carlo simulations. 
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Figure 4: Current estimates above the critical point p > p c and Jl,n > 4>a for 
L = 1000, 7 = 0.5 and 6 = 4. The fluid approximation (f^5|) and the current 
matching (|29p agree well with canonical numerics and the metastable branches 
from Monte Carlo simulations (cf. Figs. [1] and [3]) . 



4.2 Thermodynamic Limit Rate Function 

In the thermodynamic limit we expect the existence of a rate function I p that 
describes the asymptotic probability of observing a background density pbg €E 
(0, p). Precisely, we will consider the limit 

777 

N, L, m — > oo such that N/L — > p and — — ► p — pbg , (30) 
and expect that 

Tt L N [M L = m] ~ e - i - r '' ( ^) . (31) 
To derive this we define the following function for finite systems 

lL, N (m) = -- logTTi.iv [M L = m] , (32) 
which can be written as 

(L - 1) / . . N-m \ 
lL,N{m) = — (^logz ro (0) - — - log j + 

- i log W (m) + i log Z(i, iV) - ^ log L+ 

- Ilog^n 1 [^L-i=iV- 777]+0(e-^ ml " T ) . (33) 

This is derived in Appendix I A . 1 1 and is valid for any </> £ (0, oo). The first line 
in Eq. Q33p resembles the thermodynamic entropy of the background (cf. Eq. 
(|18|)). the first term on the second line is the contribution due to the maximum 
occupied site. The second term on the second line is the canonical normalisation 
and is independent of m. All the other terms will vanish in the limit (|30l) and 
we get the following result. 
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Figure 5: Thermodynamic limit rate function for a fixed system density p = 0.75 
with 7 = 0.5 and 6 = 4. I p (ph s ) is shown by a solid black line. The dashed 
lines show lL,N{m) against pbg = (N — m)/(L — 1) with N/L = 0.75, for several 
values of L calculated by exact numerics using Eq. (|2"2"j) . Note convergence to 
Ip(Pbg) (solid black) is slow and initially non-monotonic above p c . 



Theorem 1. In the limit LSU\) Ih,N{pbg) converges for all p b g G [0,p) as 

Ip(Pbg) ■= hm I L ,N{m) = s(p) - s(p bg ) , (34) 

L— s-oo 

where s is the thermodynamic entropy given by Eq. j!8\) . So for p > p c (super- 
critical case), Ip(pb g ) = for each p\, g > p c (see Fig. [5|). 



Proof. The proof follows directly from previous results on the equivalence of 
ensembles, for details see [21]. It has been shown that 

jlogZ(L,N)^s(p). (35) 

We use Eq. and choose <j> = $(pb g ) as defined by Eq. flH]). The cutoff 
to ~ N — Pbg(L — 1) diverges linearly in L for all pbg € [0, p). It follows that for 
< <f> c we have z m (</>) — > z(0), and therefore with Eq. (fig)) 

AT — to 

logz m ($(p bg )) - — - log$(p bg ) -> ^(pbg) (36) 

It follows from the asymptotic behaviour of the single site weights Eq. ([TT]) that 
-^logw(m) — ► 0, and the terms -^logi and 0(e~ TT ~y m 1 ) also vanish in the 
limit Eq. (f30|). It remains to show that logz^"* j m [Sl-i = N — to] -> 0. In 
the case pbg _ Pc we have <I>(pb g ) _ </>c and the first and second moments of 
u $(Pb ),m converge to those of v<&( Ph )■ The local limit theorem for triangular 
arrays (Thm. 1.2 in |llj ) covers the sum, Sl, of independent random variables 
whose distribution depends on the number of terms L, which is the case here 
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Figure 6: Current matching on a finite system with L = 1000 and p = 0.7 
close to ptrans(lOOO) for 7 = 0.5, 6 = 4. Left: The current in the background 
^m(pbg) and current out of the maximum g(m) plotted against pbg = (N — 
m)/(L— 1). Right: Exact numerics of Il.n using (|22[) shown in solid black and 
approximation (|38p in dashed green. Vertical lines show the correspondence 
between current matching and critical points of Il.n- The first local minimum 
corresponds to the lower stable branch of the current matching curve in Fig. 0] 
and the local maximum corresponds to the upper unstable branch. 



via the cut-off to. Convergence of the first two moments then implies 

vk7 l ^ [Sl-i = N - m] ~ , (37) 

and the result follows immediately. The case pbg > Pc can be reduced to the pre- 
vious case following [2Tj, by arranging the excess mass N — p c (L—\) in the back- 
ground among a finite number of sites so that each site contains at most to par- 
ticles. This provides a sub-exponential lower bound for [Sl-i = N — to] 
which completes the proof. □ □ 

Previous results in the thermodynamic limit imply that for N / L — > p > p c 
the background density converges to p c (see Section [2751 Eg. (fTTJ)) '). Naturally, 
the thermodynamic limit result gives no indication of the sharp transition from 
fluid to condensed putative phases or the metastable switching observed on large 
finite systems. For p > p c the system will appear to be condensed for sufficiently 
large L, since ptrans(L) \ p c as L — > oo (see Fig. [T]). Therefore the leading 
order finite size effects of Thm Q] do not describe the sharp transition or the 
metastability, as can be seen in Fig. [5] The apparent double well structure of 
Il.n is less pronounced for higher L, and the location of the global minimum 
shifts towards p c . For pbg = P, i-c. m/L — > 0, the limit of Il,n actually depends 
on the precise scaling of m, which we do not discuss here (see [4] for more 
details). To capture the sharp transition and the metastability we will replace 
the limit Eq. (|30[) by an appropriate scaling limit in section [5] 

4.3 Heuristics on Metastability 

We can understand the metastability on a heuristic level in terms of a simple 
current matching argument. lL,N(m) can be calculated efficiently for system 
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sizes L < 4000 using exact numerics (see Eq. (|2"2")l ). In Fig. [5] we see a clearly 
defined double well structure close to Ptrans(i) which accounts for the observed 
switching behaviour. We may choose <f) = & m (pb s ) instead of <f> = $(pb g ) in Eq. 
([55)1 so that the mean under the cut-off grand-canonical distribution is pb g - The 
gradient of Il,n{™) is given by a one step finite difference. If we assume for 
large L that ^~\^ m [Sl-i/(L - 1) x p bg ] ~ ^ 2 * L(r2 (which is valid in the 

scaling limit under certain conditions according to a local limit for triangular 
arrays (Thm. 1.2 in [2]), see Theorem [3] for details) and replace to — 1 by to, 
then a straightforward calculation shows 

Jl.jv (m- 1) - 7l,jv (to) « - (log$ m (p bg ) - logg(m)) , (38) 

where pbg = (-^ — m)/(L — 1). The first term on the right hand side follows 
in direct analogy with the thermodynamic limit for which d p s(p) = — log$( J o) 
(cf. |19]V and the second term is the finite difference of the maximum site 
contribution using Eq. ([5]). 

This result holds rigorously in scaling limit and even for relatively small 
systems (L rj 1000) Il.n is well approximated by Eq. ([3"5[) as is shown in Fig. 
[51 The approximation breaks down for background densities close to p since 
the cut off to becomes small and so the approximation to — 1 ps m is no longer 
valid. By the above argument we expect that solutions to the current matching 
equation (|29p correspond to local maxima or minima of the rate function. Whilst 
the current out of the maximum occupied site is greater than the current in 
the background we expect the background density to increase and vice versa. 
Therefore the first point that the two currents cross is locally stable and so a local 
minimum of the rate function. The next point they cross is a local maximum 
by the same argument, and there is another local minimum at a point close 
to p associated with fluid configurations, which the current matching argument 
predicts to be at the boundary. This is shown in Fig. [5] 



5 Rigorous Scaling Limit 

In this section we explore the finite size effects by examining the leading order 
behaviour at the critical scale on which metastability persists. As L — > oo this 
scale is given by 

N = Pc L + SpL 1 -" + o(L 1 - Q ) 

to. = (dp - Sp^L 1 " + o(L 1 ~ Q ) for p bg G [0, Sp) (39) 

for suitable a e (0, 1). We may equivalently express it as 

N/L = p c + 5pL- a + o(L- a ) 
N — to 

— — = Pc + Sp hg L- a + o(L- a ) . 



The correct scaling exponent a to capture the canonical overshoot and metasta- 
bility can be determined hcuristically by our previous current matching argu- 
ment Eq. ([2T))) . If we assume that close to p c the fluid current Eq. (f2"5)) can 
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be approximated by the first term in the Taylor expansion around p c , then 
(^m(pbg) — 1) — L~ a . The current matching equation (j2"5)l then implies 

% L~ a ~ «?(m) - 1 ~ = 6(d> - Sp^^L^-^ , 

which leads to a = j^— • 111 order to make this argument rigorous, also to 
find the transition point and describe the metastability, we study the canonical 
distribution of the background density (or equivalently the maximally occupied 
site) in the scaling limit. 



5.1 The rate function 

We examine the asymptotic behaviour of /i,jv(m) introduced in ([52)1 . We will 
see in Theorem[3]that there exists a unique 1 > j3 > such that lim L x ~^Il iv(m) 

L—>oo 

is finite and non-zero at the critical scale, we define 

4?(<Vb g ) = = I™ L l -Pl L , N {m) , (40) 

where N and m are given by Eq. ([M]). Since I p (pbg) = for pbg > Pc (see Sec. 
4.2) this definition implies 

ttl,n [Ml = m] ~ e ~ Lf>I ^p {5phs) (41) 

as L — > oo, we recall that if Ml = m the background density is N L ~^ and 5p 
and (5/Obg are given by Eq. ()39l) . We begin by describing the scaling limit for the 
current in the fluid bulk and out of the maximally occupied site. 

Theorem 2 (Current scaling limits). For the scaling limit with a = y-p^. 

the current out of the most occupied site and the average current in the fluid 
background are asymptotically given by 



provided -^rSptg < tj^ 



!, 



(l-y)(6p-6p bg )-y ■ 

Proof. The result for the current out of the condensate g(m) is immediate since 
to = (Sp - 5p hg )L 1 " a + o^ 1 ""), g{m) = 1 + bm^ and (1 - a)j = a. 

The proof of the fluid result follows from a Taylor expansion of the trun- 
cated density function R m (<j)) introduced in (|26[) . details on the expansion are 
contained in Lemmas [T] and [2] in the appendix. For convenience we write <f> = e M , 
then 



for all p, £ [0, , 1 _ \(g b _g — w )• The result follows directly since 



R m (e» L ~ a )=Pc + <J 2 c pL~ a + o{L- a ) 

b • 

e' ii " a =* m (i2 m (e' u '"")J 
=> 1 + pL~ a + o(L- a ) = <P m (p c + a 2 c pL~ a + o(L~ a )) 
and choosing p so that a 2 p — <5pbg- d d 
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Although the proofs in this section are restricted to -4-<5phe < r\ — vrr — 5 — v 

o 1 ai (l-7)(5p-(5pbg) 

this includes the point of intersection of the two currents (cf. Fig. [S]) and we 
discuss a possible extension to the whole region 5p\, g G (0, Sp) at the end of this 
Subsection. Theorem [2] states that in the scaling limit the average current in 
the fluid background and the jump rate out of the condensate converge to the 
critical current <f> c = 1 on the same scale as the density N/L converges to p c . 
So for N and m at the critical scale, given by Eq. ((39]) . (1 — ^ > m (pbg)) L a and 
(1 — g(m)) L a converge to unique functions of the scaled variables Sp and Sphg- 
We now show that in the scaling limit the rate function is exactly the integral 
of the difference in the two limiting currents (up to normalisation). 



Theorem 3 (Scaling limit rate function). In the scaling limit \S9\) . with a 
j-pj, the rescaled rate function J^0| ) with j3 = -j-^ converges as 



ifftfptg) = lim L 1 ~' 3 / L , w (m) 



L 

,2 



2al 1-7' y/ re(o,<5p) |_ 2<7<f I-7' 

provided ±Sp bg < {1 _ lKS b p _ Spbgr , ■ 
Proof. Firstly consider the unnormaliscd measure 



Hl,n = — L , 1N ^L,N , (42) 



Z(L,N) 
^(1) 

for which it follows from definitions in Section 12.21 and Eq. (|32|) 

ilog/ii.jv [M L = m] = -L 1 ^ : I L , N (m) + -^logZ(L,N) - L^log^l) 



(43) 



Re-writing Il,n using Eq. (j3"3")l we find 

log/ii,jv [M L = m] = - ^jf^ (logz m ($ m (p£ 8 )) - p^log* m (p^ g )) + 
+ L 1 "' 3 logz(l) - log to(m) - log£+ 

- 77 lo g<;U s ),™ [5l - i=7V ^ m] ~ o(1) ' (44) 

where we used the shorthand p^ g = (N — m)/(L — 1). By applying Lemma 
[51 given in the appendix, the variance under the cut-off distribution v$ m ^ p L j m 

converges to the variance of the critical measure for -^Spbg < (i-~f)(Sp-8p b p ■ 
By the local limit theorem for triangular arrays (Theorem 1.2 in [11] ) we have 

l v = -/V - ml ; 1 (45) 



since — > p c and E„^ L [n x ] — > p c as L — > 00 with the scaling 

So the contribution of this term to (|44p vanishes. Also since /3 > we have 
i - ' 3 log L — >• as L -> 00. 
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The condensate contribution to Eq. (|4"4"f is given by 

(Sp - Sp^) 1 ^ (46) 



1-7 



where we have used the asymptotic behaviour of the single site weights (Eq. 
GU) and (l-a)(l-7) = 0. 

By applying Theorem [5] to ^m(Pb g ) an d taking Taylor expansion of z m ac- 
cording to Lemma [5] (appendix), 

L 1 -? (p c + ^ bg L-«) log $ m ( Pc + <5pb g L- a ) = %^> c + ^ + o(l) (47) 

O C 

and 

log*™ (*„,(/&)) = (48) 
= L 1 -? log ^z(l) + ^iiTV(l) + ^L- 2 «z"(l) + o(iT 2 ^ 

= log z(l) + %L> C + ^| + (1). (49) 

2cr 2 

since (1 — 0) = 2a. Combining Equations (|45D to (|48[l in (l4"4l implies, 

-i log/ii.jv [M L = m] = + -^-(fy - ^ bg )^ 7 + o(l) (50) 

in the scaling limit as L -> oo, provided ^5pbg < (i_ 7 )(^p-g Pb w • 
By definition 

TV 

Z(L, N) = z L (l) J2 Vlm [M L = m] . 

m=l 

So we may bound Z(L, N) as follows, 

z L (l) max{p L .N [M L = to]} < Z(L, N) < Nz L (l) ma,x{p L ,N [M L = m]}. 

Applying Eq. ([50)1 the asymptotic behaviour of L~P log Z(L, N) in the scaling 
limit is given by 

L-P log Z(L,N) =L 1 -?lo g z(l)- + -^(^-^1+0(1) , 

re(o,8 P ) I la* 1 — 7 J 

(51) 

provided that the global minimum of — jj log pl, n [Ml = to] is attained for 
-^zSpbg < ^ 1 _ 7 ^gp_g Pb ) 7 in the scaling limit. This follows from results in [3]. 
Combining Equations ([50)1 and ([5Tj) the result follows from definitions Eq. ([55)1 
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and Eq. @2J 



- inf {^— + -^—(6p-r) 1 -A+o(l) . 

□ □ 

The result may break down for ^<5/9bg > (i_ 7 )(gp_^ Pb re due to the large 
probability of a second sub-condensate forming on the same scale as the con- 
densate. The conditions under which this occurs could be found by conditioning 
on a maximum site occupation and using the same methods as for Theorem _J| 
to find the distribution over the second highest occupied site. By using this 
method and iterating when necessary, the results here could be extended to 
the entire region bp^ g € (0,<5/o). However our results already cover the relevant 
critical points of the rate function, the fluid minimum, condensed minimum and 
local maximum that constitutes a potential barrier. 

5.2 Current Matching and Overshoot 

Theorem [3] together with Theorem [2] imply that in the scaling limit given by 

_____ /Q 1 . 

([55)1 the derivative of the scaling rate function Ig ' is given by the difference in 
the current out of the condensate and the average current in the background. 
This is made precise in the following result, which is a rigorous version of the 
current matching argument of Section [4] 

Corollary 1 (Current matching). 

^<>. 3 ) = ^-^_^ (52, 

and the two terms on the right-hand side are exactly the limiting current curves 
from Theorem above the critical point. 

This implies that 

d 6pbJsp ( s Pbg) -> -°o as (fog -> , (53) 

(2) 

so that Ig always exhibits a boundary minimum at bp^g = Sp. 

Depending on the value of Sp the rescaled rate function has one of three 
qualitative forms characterised by the number of its extreme points, or cquiv- 
alently the number of roots of (|52[) , and the position of the global minimum. 
With the threshold 

coh,b) = — — (cr c 7&) , (54) 

Eq. (|52[) has no real roots for bp < Co, exactly one for bp = Co, and two for 
bp > Co- The latter correspond to a local minimum at i5pbg = tq € (0,bp) 
and a local maximum at <5pbg = f*, with ro < r* (cf. Fig. [7] on the right). As 

(2) 

usual, minima of Ig correspond to metastablc phases, and the depth of the 
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local minimum at tq as compared to the one at the boundary determines which 
of the phases is stable (i.e. corresponds to the global minimum). So cq marks 
the threshold above which (for Sp > cq) there exists a metastable condensed 
phase. Stability of the phases changes when both minima have the same depth, 
the density where this is the case is given by ptrans{L) introduced in Section 3. 
In the scaling limit (|3l?l) this behaves as 

Ptrans(L) = p c + 5p tr an S L~ a + o(L~ a ) , 

and from Theorem [3] and Corollary [T] we get the explicit expression 

i 

Note that both Co and Sptrans are increasing with 7 and b. Their ratio simplifies 
to 

1/ 2 \ 1/(1+7) 

<Wan S /C0 = 2\'l^~ ) ' ( 56 ' 

which also increases monotonically with 7. In particular, Sptrans > Co for all 
7 > which implies existence of an extended metastability region. 

We can now summarise the critical behaviour of the system at scale (f3U| in 
three cases, which are also illustrated in Fig. [7] (right). 

Case 1 (Sp < Co). There exists a unique minimum of the scaling rate function 
at Sphg = Sp, 

4?(M = o, 

corresponding to a stable fluid phase. The most likely background density under 
the canonical measure is S b L e /(L - 1) = p c + SpL a + o(L- a ). 

Case 2 (co < Sp < 5/Otrans)- There are two local minima of the scaling rate 
function 

corresponding to a stable fluid phase with <5pbg = Sp and a metastable condensed 
phase with <5/9bg = ?*o- The most likely background density under the canonical 
measure is still S h L s /(L - 1) = p c + 8pL~ a + o{L~ a ). 

Case 3 (Sp > Sp tla , ns ). There are two local minima of the scaling rate function 

/g ) (M>4?M = o, 

corresponding to a metastable fluid phase with Spb g = Sp and a stable condensed 
phase with Sph s = ro- The most likely background density under the canonical 
measure is now S^ s / (L — 1) = p c + r a L~ a + o(L~ a ) and a finite fraction of the 
excess mass condenses on a single lattice site with m = (Sp — ro)L 1 ~ a +o(L 1 ~ a ). 

The value of Sptrans, and for Sp > (Jptrans the position of the global minimum, 
agree with recent results on the distribution of the maximum at the critical 
scale [3]. The expected value of the maximum site occupation (equivalently 
background density) under the canonical measure can be interpreted in terms 
of the stable solution of the current matching argument as described above. 
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Figure 7: Scaling limit for 7 = 0.5 and 6 = 4. Left: The thick black line 
shows the scaling limit result including metastable branches (dashed). Thin 
coloured lines show the rescaled canonical currents from exact numerics, slowly 

approaching the scaling limit as L increases. The transition point is given 

(2) 

by <5ptrans- Right: I Sp for various values of S p. Solid black corresponds to 

Sp = <5ptrans where the depth of the two local minima are equal. The green 

(2) 

dashed-dot-dot curve shows 1$ ' for op < Co, blue dot-dashed curve corresponds 
to cq < 5p < (Sptrans and red dashed Sp > S p trans- 



Corollary 2 (Canonical current overshoot). Under the conditions of Theorem 

H for Sp < Sptrans 



3L,N-l = L- a %+o{l) (57) 



and for Sp > 5p tra ns 



bL~ a 

.]L,N - 1 = T , rj + o(l). (58) 

{op- r ) 



Proof. Follows directly by expressing Jl,n in terms of the ratio of canonical 
partition functions Eq. pH|) and applying Eq. ([51]) . □ □ 

This result is illustrated in Fig. [7] (left). At the transition point Sp = (5/Otrans 
we must examine higher order terms to demonstrate that the system will asymp- 
totically be in a condensed phase, see [3] and section [B] 



6 Estimating the Lifetime of Metastable Phases 

Having established the existence of metastable phases we can use our previous 
results to get an estimate of their lifetime. Previous studies include results on 
the dynamics of fluctuations in a finite system |22) and the dynamics of the 
condensate [TTJ [TBI IS] • O ur approach here follows mostly the one in [5D] where 
metastability in a different zero-range process has been studied, and |18j where 
a random walk argument was used to find the characteristic time of the motion 
of the condensate. 

In this section we consider the dynamics in the scaling limit defined by Eq. 
(pii?)) . In Section 5.2 we have seen that for Sp > cq the system exhibits two 
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metastable phases, we assume in the following that we are in this region. So 
for sufficiently large L we know -Zl,at(to) exhibits a local maximum between the 
fluid and condensed minima. Let to* be the local maximum of Il.n which is 
well approximated by the closest integer to the largest root of Eq. (j2"9")l (see 
Section [4]). In the scaling limit N ^"f = p c + r*L~ a (see Section [5]). We define 
the metastable fluid (condensed) phase as all states for which the maximum site 
occupation is below (above) m* , i.e. 



^COIIC 



cond ~ n L , N [-\M L (r,) > m*] . (59) 



Since the process spends very little time close to the maximum of Il.n, the 
precise choice of to* is not relevant for results in the scaling limit (cf. also [7] 
for a slightly different approach). The lifetimes of the two metastable phases 
can be expressed in terms of the following hitting time, 



We take 



r L :=mf{t>Q:M L (r)(t)) = m*} 



Tfluid := 2E 7r fiuid e£t [tl] and T cond := 2E 7r co„d e£t [t l ] , (60) 



where the expectation is with respect to the dynamics Eq. @ with initial dis- 
tribution given by irfj 11 ^ or 7r™ £ d . The factor of 2 comes from the fact that once 
to* is reached the process can still return to the same metastable phase with 
probability 1/2. 

By ergodicity of the process (rj(t) : t > 0) on SIl,n the ratio of lifetimes is 
directly related to the stationary distribution, 



T flui d _ TT LiN [M L (ri) < to*] £™ =0 e -LiL,N(k) 



T C o„d itlA m l(v) > m*} Y.k= m *+i e ~ LlL - N(k) 
Applying the scaling limit result of Theorem [31 we get for large L 

T nuid f 5 " - L ^< 2 > 



(61) 



2~cond 



e L 1& p ^ dx e L 1§ p ^ dx J 

V2^L 



(62) 



where the second line follows by a saddle point approximation of the denom- 
inator at the local minimum ro that corresponds to condensed configurations, 
and by expanding II 1 (x) in the numerator at the boundary minimum (corre- 
sponding to fluid configurations), keeping the leading order singular term. To 
leading order in the exponent this result implies that the ratio of the lifetimes 
is given by the relative depth of the minima, as expected. The lower order 
term \f~L in the denominator implies that at the transition density <5ptrans in the 
scaling limit the condensed phase is stable and the fluid phase is metastable, in 
accordance with results in [3J. However the scaling with L is difficult to verify 
in Monte Carlo simulations, since system sizes have to be relatively small (at 
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Table 1: Lifetime of metastablc phases for L = 1000 and N = 695, with 7 = 0.5 
and 6 = 4. Using approximation (|63p compared to measurements from Monte 
Carlo simulations on a fully connected lattice (MF), a one dimensional lattice 
with p = 1 and with p = 3/4. 



Approx. ((631) 



MF 



lDp=l lDp = 3/4 



Jfluid 
-*cond 



1.61 x 10 5 
1.51 x 10 5 



1.65 x 10 5 
1.52 x 10 5 



1.21 x 10 5 
1.07 x 10 5 



2.05 x 10 5 
1.94 x 10 5 



most of the order of L = 1000) to get good statistics on switching times. In this 
regime higher order finite size effects still play a role and can affect the location 
of Ptrans(i) and therefore the relative depth of the minima of Il,n- 

We estimate the two lifetimes by approximating the number of particles in 
the condensate by a continuous time random walk. Since it is difficult to validate 
the scaling with L, as discussed above, we will demonstrate the validity of our 
approach by estimating the lifetimes directly for finite systems. The actual 
lifetimes are not only related to the potential barrier given by the maximum 
of Il.n, but also depend on the underlying dynamics. The occupation of the 
maximum (Ml (77(4)) : t > 0) is a non-Markovian, stationary, ergodic process 
with state space f2* = {\N/L], \N/L]+1, . . . , m*, m*+l, . . . , N} and stationary 
distribution tt*(to) = e --"x.,JvM. Since the process exhibits only single steps, 
it can be approximated by a continuous time random walk on fi* (the validity 
of this Markovian assumption is discussed later) where a particle leaves the 
maximum (to — > to — 1) with rate g (to). The rates corresponding to to — ¥ 
to, + 1 are then fixed by the stationary distribution 7r* . In the scaling limit this 
implies that particles enter the condensate with rate given by the current in 
the background, see Corollary [TJ The lifetime of the metastable fluid phase is 
approximated by the mean first passage time of the random walk starting in 
state \N/L~\ to reach to*. For the metastable condensed phase it is given by 
the mean first passage time starting at N, and both have a factor of 2 in front 
due to the possibility of reaching to* and not actually switching phase. These 
can be calculated using standard techniques (see e.g. [35] for the discrete time 
analog), 

m* x / . ^\ m* ^ z—2 

Ifluid ~ 2 V + 2 V V 7T*(j) 

q[i)ir*(i) t—^ q(i)iT*(i) ' 



N N , N 



- cond 



2 E 4r + 2 E TTTm E ^0')- (63) 

i=m*+l 3V ; i=m* + l lVW w j=i-l 



These two approximations can be calculated numerically for reasonably large 
systems (up to L ps 4000) since we can calculate 7t*(to) = e~ L/ i.«( m ) exactly 
using Eq. (|2"2")l . Applying the scaling limit result and using saddle point approx- 
imations we recover an Arrhenius estimate with a prefactor for the lifetime of 



22 



the two metastable phases, 



1 



[ e - Lf>I %Mdydx 

J X 



T, 




cond 

Jo 

„ t£ e ^(4?(r*)-I$(r ))_ (64) 

Here r* is the location of the local maximum of the rate function, r$ is the local 
minimum corresponding to condensed configurations and Sp is the boundary 
minimum corresponding to fluid configurations. To leading order the lifetime 
of the metastable fluid (condensed) phase is given by the height of the local 
maximum of the rate function (potential barrier) above the fluid (condensed) 
minimum. This result agrees with the ratio of the lifetimes already discussed. 
However as with the previous results the asymptotic form is difficult to validate 
due to higher order finite size effects, the exact expression Eq. ((63)) is therefore 
more appropriate for relevant system sizes. 

We expect the above random walk approximation to be accurate if particles 
exiting the maximum equilibrate in the bulk before returning to the condensate. 
This condition is best fulfilled for a fully connected lattice (mean-field geometry) 
and we expect it to be a reasonably good approximation for one dimensional 
totally asymmetric systems since particles have to pass through the whole fluid 
bulk before returning. For partial asymmetry and in higher dimensions a re- 
turn without penetrating the fluid bulk is possible. To a first approximation 
the lifetimes have to be multiplied by the inverse probability of the event that 
particles escape into the bulk, since only such particles have a chance to equili- 
brate and contribute on the right scale. This pre-factor can be estimated using 
another random walk argument. If a particle jumps to the right with probabil- 
ity p and to the left with probability q where p + q = 1 and p ^ 1/2 then the 
probability of a particle reaching some macroscopic distance before returning 
to the condensate is asymptotically \p — q\. We therefore expect the lifetime to 
increase by a factor of l/\p — q\ for partial asymmetry. As a special case, for 
symmetric systems in one dimension this leads to an increase of lifetimes by a 
factor of L by the same argument. For a detailed investigation and validation 
of this Markovian ansatz see [IS]. We find that this argument gives a good 
approximation for fully connected lattices, and totally asymmetric jumps in one 
dimension, see Table [TJ We observe that the one dimensional totally asymmet- 
ric case is a factor of approximately 1.4 faster than the mean- field case, which 
is due to internal structure in the fluid background, and has also been observed 
and discussed in [IS]. For partial asymmetry in one dimension with p = 3/4 the 
process is slower than the totally asymmetric case by a factor of approximately 
1.8 supporting the arguments above which predict an increase by a factor of 2. 

The characteristic time for the motion of the condensate in the thermody- 
namic limit, N/L — > p > p c , can be approximated by considering Eq. (|39[) with 
a sequence of Sp increasing like L a . The corresponding rcscalcd location of the 
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condensed minimum r<j — > 0, since the condensed phase has limiting background 
density p c . Our rigorous results do not technically hold in this limit, however 
our estimates of the lifetimes of the metastable phases agree with the character- 
istic times found by Godreche and Luck [TBI HZ] • Our results demonstrate that 
close to the transition point in finite systems (i.e. Sp w 5/Otrans) the condensate 
typically moves via the system entering the fluid phase (cf. Fig. [3]), and therefore 
its new position is expected to be chosen uniformly at random on A^. For Sp 
large, however, the lifetime of the fluid phase becomes small and the condensate 
can re-locate whilst the system remains in the condensed phase, via the mech- 
anism described in [TS] . A sub-condensate starts to grow in the fluid bulk, and 
the potential barrier the process has to cross for condensate motion is associ- 
ated with the probability of the excess mass shared equally between two sites. 
This is the relevant mechanism in the thermodynamic limit. For finite systems 
one can estimate which of these re-location processes dominates, by considering 
the distribution of the second largest site using the same techniques as in the 
previous section. It has recently been shown that the second mechanism via a 
sub-condensate can lead to a non- uniform relocation of the condensate [BJ. 

7 Discussion 

We have shown that a prototypical class of zero-range processes, that are known 
to undergo a condensation transition, exhibit large finite size effects including 
a metastable switching phenomenon for a large range of system parameters. 
We have characterized this behaviour using the background density as an order 
parameter, by rigorously deriving a large deviation rate function for its distribu- 
tion in an appropriate scaling limit, which shows a double well structure. These 
results agree with recent work on the zero-range process at the critical scale 
[3J, which we extend by establishing metastable fluid and condensed phases. 
Our methods give rise to a simple interpretation of these results in terms of a 
stationary current balance between the fluid background and the condensate. 
All results presented here, except those in Section [BJ concern properties of the 
canonical stationary distribution and are therefore independent of the geome- 
try or dimension of the lattice, so long as it permits homogeneous stationary 
distributions. 

Our results allow us to estimate the lifetime of the two metastable phases 
using a heuristic random walk argument. The estimates agree with previous 
studies of the dynamics of the condensate in the thermodynamic limit |18j . We 
show that in the critical scaling limit the lifetime of the two phases is growing 
with system size (as a stretched exponential) and derive the appropriate Arrhe- 
nius law including prefactor. We have also demonstrated that on finite systems 
the re-location dynamics of the condensate is dominated by switching to the 
metastable fluid phase for a large parameter range, rather than the mechanism 
via a second sub-condensate assumed in previous results, which is relevant in 
the thermodynamic limit. 

The zero-range process provides an interesting example where the thermody- 
namic limit fails to capture essential parts of the dynamics, and contradicts the 
usual expectation that finite systems should behave in a smoother fashion than 
the limiting prediction. Although the aim of this paper is not to discuss partic- 
ular applications in detail, the finite size effects described here are particularly 
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Figure 8: The current overshoot in the scaling limit for 7=1 and 6 = 4. 
The thick black line shows the scaling limit result and the rescaled canonical 
currents are shown for various system sizes by thin coloured lines, they slowly 
approaching the scaling limit as L increases. The transition point is given by 
Co (7, 6) which is given in [3]. 

important for a proper understanding of various real condensation phenomena. 
To name two relevant examples, the condensing zero-range process has recently 
been applied as a simplified traffic model as well as to describe clustering phe- 
nomena in granular media [351 13H1 133 130] • Both of these applications exhibit 
metastability phenomena, and typical system sizes are of order 10 2 — 10 3 , the 
region where metastability effects in our analysis are most relevant. Since our 
results hold for a large generic class of zero-range models, they suggest that the 
apparent metastability in many applications is not necessarily the result of a 
particular choice of the jump rates, but rather a generic finite size phenomenon. 
While the results in the scaling limit hold for all zero-range processes with the 
same tail behaviour as the rates (TTJ), the strength of the effect for relevant system 
sizes will of course depend on the details of the particular application. Traffic 
modelling is an example where metastability is particularly pronounced [41] . 
and it is an interesting question to investigate in detail the applicability of the 
zero-range process and our results in that area. 

Using the methods discussed in this paper, a direct application of the findings 
in [3] allows us to extend our results on the current overshoot to zero-range 
processes with rates ([T]) where 7 = 1 and b > 3, see Fig. [5J The critical scale 
in this case turns out to be of order y/L/ logL, and the system also shows a 
discontinuous current overshoot in the scaling limit. However, as can be seen 
in Fig. [5J convergence is very slow and there are no metastable phases on the 
critical scale. Although metastability might occur on higher order scales, it will 
hardly be relevant in any finite size simulation or application. 

The scaling limit rate function in Theorem [3] could also be derived using 
a slightly modified version of the saddle point methods applied in [Ml [15l [32] . 
However, the direct approximation of the fluid and condensed metastable phases 
used in this paper seems more intuitive for our purpose, and can be made 
rigorous with less effort than a saddle point computation. 



— Scaling limit 

— L=10000 
L = 2000 

— L = 1000 
L = 500 
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Appendix 

A.l A formula for Il,n 

In Eq. (j52"|) we introduce Il,n{iti) = — j^logiTL,N [Ml = m] for the canonical 
probability that the maximum site occupation is m. The following upper bound 
over counts configurations in which more than one site contains m particles 

ttl,w [M L = m] < ' N > , (65) 

where X = {rj : 771, ... , 771,-1 < m, YLx=x Vx = N — m}. The following lower 
bound does not count any configurations in which more than one site contains 
m particles, 

ttl,n [M L = to] > ~Z{LN) ' *• ' 

where X = {77 : 771, . . . , ?7l-i < to — 1, 2x=i J te = N — m}. It follows that 

«l,n [M L = m] = zIVn] °^ Ler ~ m )) ( 67 ) 

as L,to — > 00. We may write the last term in the numerator in terms of the 
cut-off grand canonical measure and get 



kl,n [M L = m] = 

^\i-(A r -m)..I-l[P j\r _1 

(l-0(Le-^ ml_T )) , (68) 



Lw{m)z^4>)r {N - m) ^[SL-x=N-m] , __L m i-T, 



Z(L,N) 

which holds for all <fi £ (0, 00) (cf. Section 4.1). Finally, taking logarithm, 
-v log n L .N [M L = m] = - ^ - ( log z m (<f>) - N _™ log <pj + 

- - log w(m) + - log Z(L,N)-- log i+ 

- \ log^" 1 [5£_i=JV - m] +0{e-^™ 1 -'). (69) 

A. 2 Convergence of Moments under the Cut-off Ensemble 

Lemma 1. Consider two sequences m n G N, fj, n G (0,oo), such that m„ — > 00, 
Hn — > and tn^fin — > C as n — > 00 with C <G [0, j^zr)- Then for each i € No 

limsup ^2 k 1 w(k)e ,lnk < 00. 
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Proof. From the asymptotic behaviour of w(k) Eq. (fTTj) we know there exists a 
C 6 (0, oo) such that ui(fc) < C e~ fcl 7 for all k. So, 

^ fc I w(fc)e' I " fc < C k l e^ kl ~ 1+ ^ k . (70) 

k<rn n k<.m n 

Fix < e < (j^- - C). Then there exists n e N such that /z„ < (C + e)m r 7 7 for 
n > n. It follows that \i n < (C + e)fc~ 7 for all k < m n and n > h. By applying 
this upper bound on fx n we can bound the sum above uniformly in n > n as 
follows, 



k l w(k)e^ k < C J2 fc ' e 



-' + (C+e)k 1 



fc<mn k<m n 

/■oo 

< / x * e T^^HC+e)^ < QO _ (71) 

Jo 

□ □ 

Lemma 2. Under the assumptions of Lemma{J\we may bound the remainder in 
the Taylor expansion of each moment of the cut-off measure ^ >jmn introduced 
in i23\) . Writing <j) n = e^ n we have 

z m n {^ n ) = + z'{l)n n + o(fin), 

Rm n (e Mn ) = p c + 0-\ll n + o(/i„), 
and . m „ H] = E tff , B fol] + Und^^ [Hi] + o(Mn). 

Proof. For each n £ N the i moment of the cut-off ensemble v^> n ,m„ is given 
by a finite sum. We let, 

/m n (M)= E *Mfc)e" fc - (72) 
then /„ ln £ C°°(M+,R+). By Taylors theorem expanding around /x = 0, 

fm n (fhi) = /m„(0) + /m„(0)/i„ + h n (ji n ), 

where the remainder term may be expressed as 

(2) 

/m„ is non-negative and increasing and so the remainder is bounded by 

< K{n n ) < fx n . 

By Lemma [TJ we know that there exists U > such that /m„(e M ") < J7 for all 
n £ N. It follows that, 

/mje"") = /m„(l) + C B (l)A*n + 

The result then follows by considering (/^(l) - / m „(l)| - Efc> m „ k l e T ^ kl " ,+flnk 
which is o(/i rl ) since fJ/ n m% — > C, so the tail of the sum converges as a stretched 
exponential, this allows us to replace the expected values under the cut-off with 
/i = with the expected values under the true critical measure. □ □ 
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